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Abstract 

The paper deals with a theoretical investigation of the peristaltic transport of a couple 
stress fluid in a porous channel. The study is motivated towards investigating the phys- 
iological flow of blood in the micro-circulatory system, by taking account of the particle 
size effect. The velocity, pressure gradient, stream function and frictional force of blood 
are investigated, when the Reynolds number is small and the wavelength is large, by using 
appropriate analytical and numerical methods. Effects of different physical parameters 
reflecting porosity, Darcy number, couple stress parameter as well as amplitude ratio on 
velocity profiles, pumping action and frictional force, streamlines pattern and trapping 
of blood are studied with particular emphasis. The computational results are presented 
in graphical form. The results are found to be in good agreement with those of Shapiro 
et. al [52] that was carried out for a non-porous channel without consideration of couple 
stress effect. The present study puts forward an important observation that for peristaltic 
transport of a couple stress fluid during free pumping, flow reversal can be considerably 
controlled by suitably adjusting the couple stress effect of the fluid/Darcy permeability 
of the channel. It is also possible to avoid the occurrence of trapping, by reducing the 
permeability. 
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1 Introduction 

Peristalsis is a natural mechanism of pumping that is observed in the case of most physiological 
fluids. In the transport of some other fluids also peristaltic behaviour is observed. This be- 
haviour is usually associated with a progressive wave of area contraction or expansion along the 
length of the boundary of a fluid-filled distensible tube. Peristaltic transport occurs in various 
physiological activities, for example, in the flow of urine from kidney to the bladder, in the move- 
ment of food material through the digestive tract, in flow of fluids through lymphatic vessels 
as well as in semen movement in the vas differences, in the movement of bile from gall bladder 
into the duodenum and spermatozoa inside the ductus efferentes of the male reproductive tract 
and cervical canal, in flow of ovum in the fallopian tube, the movement of cilia and also in the 
flow of blood through small blood vessels. This phenomenon is also applied in the propulsion of 
some industrial fluids. 
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Nomenclature 




a 


Wave amplitudes 


j 
a 


Half-width of the channel 


Da 


Darcy number 


ti 


Vertical displacement of the wall 


k 


rermeabihty parameter 


a, 7] 


Couple stress parameters 


r> 
r 


Fluid pressure 


Q 


Flux at any axial station 


R 


Reynolds number of the fluid 


t 


1 lme 


X, Y 


Rectangular Cartesian co-ordinates 


T T T r 

u,v 


Velocity components m A,Y directions respectively 


m 1 ,m 2 


Constants defined in equation (18) 


x 




Wave number 


A ^ 

i\p 


Pressure difference between the channel ends 


e 


Porosity parameter 


A 


Wave length of the traveling wave motion of the wall 


n 


Viscosity of the fluid 


V 


Kinematic viscosity of the fluid 





Amplitude ratio 


P 


Density of the fluid 



This mechanism is also used in many biomedical appliances, such as finger pumps, heart-lung 
machine, blood pump machine, dialysis machine and also in industries for the transport of 
noxious fluid in nuclear industries, as well as in roller pumps. For this reason in current years, 
studies of peristaltic movement have been receiving growing interest of scientific researchers. 
Much of the early literature of theoretical investigations, arranged according to the geometry, 
the type of fluids, the Reynolds number, wave amplitude, wavelength and the wave shape, 
along with an account of experimental studies on peristaltic transport have been reviewed by 
Srivastava and Srivastava [1]. Some of the important theoretical studies on peristalsis have been 
discussed by Eytan et al. [2], Jimenez- Lozano et al. [3], Misra et al [4, 5, 6, 7, 8, 9, 10, 11], 
Maiti and Misra [12], Usha and Rao [13], Mishra and Rao [14], Rao and Mishra [15], Park et al 
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[16] as well as by Akbar and Nadeem[17]. A few of these studies have been carried out by using 
the lubrication theory by neglecting the fluid inertia and wall curvature without any restriction 
of the wave amplitude. Some other analyses are based on the consideration of small peristaltic 
wave amplitude, where the Reynolds number is arbitrary. Antanovskii and Ramkissoon [18] 
studied the peristaltic transport of a compressible viscous fluid through a pipe with the help of 
lubrication theory when the pressure drop changes with time, by taking into account the wall 
deformation of the pipe. 

Bergel [19] observed that the capillary walls are surrounded by flattened endothelial cell 
layers which are porous. Dash et al. [20] considered the Brinkman equation to model blood flow 
in a coronary artery in the pathological state when clogging of blood occurs in the lumen of the 
artery. They assumed the clogged region as a porous medium and assumed the permeability to 
be const ant /varying in the radial direction. Misra et al. [21, 22, 23, 24, 25, 26] studied different 
aspects of blood flow in arteries in normal/pathological states. Some models on blood rheology 
have been discussed by Chien and Skalak (cf. [28, 29]). 

It is well known that most physiological fluids including blood behave as a non-Newtonian 
fluid. It has been suggested that blood flow behaviour in small vessels (of diameter less than 
0.2 mm) at low shear rate (< 20sec~ 1 ) can be represented by a power law fluid [30, 31]. It was 
reported by Blair and Spanner [32] that blood can be represented by Casson's model for moderate 
shear rate flows. The non-Newtonian behaviour of blood is mainly due to the suspension of 
red blood cells in the plasma. When neutrally buoyant corpuscles are contained in a fluid 
and there exists a velocity gradient due to shearing stress, corpuscles have rotatory motion. 
Furthermore, it is observed that corpuscles have spin angular momentum, in addition to orbital 
angular momentum. As a result, the symmetry of stress tensor is lost in the fluid motion that 
is subjected to spin angular momentum. The fluid that has neutrally buoyant corpuscles, when 
observed macroscopically, exhibits non-Newtonian behaviour, and its constitutive equation is 
expressed by a stress tensor [41]. The radius of gyration of the corpuscles for such fluids is 
different from that of the fluid particles. Their difference produces couple stress in the fluid. It 
is known that such a fluid has spin angular momentum in addition to the couple stress effect. 
The importance of consideration of couple stress effects in studies of physiological and some 
other fluids was indicated by Cowin [33] and Beg et al [34, 35]. 

Studies pertaining to the couple stress fluid behaviour are very useful, because such stud- 
ies bear the potential to better explain the behaviour of rheologically complex fluids, such as 
liquid crystals, polymeric suspensions that has have long-chain molecules, lubrication as well 
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as human/sub-human blood [41, 36, 37, 38]. Valanis and Sun [39] as well as Popel et al.[40] 
remarked that couple stress fluids constitute a special class of non-Newtonian fluids that take 
account of the particle size. The couple stress theory presented by Stokes [41] is quite suitable to 
study blood flow in micro- vessels by taking the size of the erythrocytes into account. Shehawey 
and Mekheimer [42] analysed the flow of a couple stress fluid for any arbitrary Reynolds number 
and wave number, considering the wave amplitude to be small. 

It may be mentioned that consideration of porosity is very much necessary to properly 
explain the fluid dynamical process that occurs in different parts of the body, such as vascular 
beds, lungs, kidneys and tumorous vessels. The importance of consideration of Darcian drag 
effects in studies related to flow of blood and some other fluids through porous media has been 
discussed by different researchers [43, 44, 45, 46, 47, 48, 49, 50]. Moreover, in many biomechanical 
studies, porosity of the media has significant influence on the transport of fluids. This applies 
more particularly to vessels impeded by clots and also to highly perfused skeletal tissues, tumors 
and to soft connecting tissues [51]. Some further discussion on vessel porosity and peristaltic 
flows of physiological fluids is available in several scientific publications [52, 53, 54, 55, 56, 57, 58]. 

Being motivated by the observations reported in the above mentioned studies, we have 
undertaken here a study that concerns peristaltic flow of blood in a porous bed. With an aim 
to take account of the particle size effect, blood has been modelled as a couple-stress fluid. The 
problem has been analysed by using lubrication theory [52]. The analysis and the results of 
the study are particularly applicable to the peristaltic transport of blood in coronary arteries of 
smaller dimensions in the pathological state, when the lumen of the segment of a small blood 
vessel turns into a porous medium due to the presence of numerous blood clots, or when arterial 
clogging takes place by deposition of fatty plaques of cholesterol in the arterial lumen, and also 
in cases where numerous tumors are grown inside the lumen due to excessive cell division. 

The results of the present study will serve as a reasonably good estimate of various fluid 
mechanical parameters for peristaltic transport of blood in small blood vessels in a pathological 
state. Since flow behaviour in an axi-symmetric vessel resembles that in a channel, we have 
studied here two-dimensional channel flow of the fluid. The results of this study will be applicable 
to blood vessels in the micro-circulatory system without any restriction. 

The study bears promise of multi-fold important applications. In the realm of physio- 
logical fluid dynamics, in many situations it is required to have estimates of a variety of 
fluid mechanical variables when some physiological fluid has to pass through porous struc- 
tures, particularly in pathological states. The study will also have an important bearing in 
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examining the flow in a vessel when the luminal surface of endothelial layer is attached with 
glycocalyx, which contains a series of micro-molecules and adsorbed plasma proteins. The 
present study can also find important applications in artificial bio-processors that have syn- 
thetic porous surfaces. Since blood flowing in small vessels has got relevance with flow in 
porous media, the study is applicable to blood flow in micro-vessels more particularly in sit- 
uations where couple stress effects of blood are prominent. In addition to all these, the anal- 
ysis presented here should find applications in polymer industries, where it is often required 
to examine the flow behaviour of different types of fluids passing through porous structures. 
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Fig. 1: A physical sketch of the problem 



2 Mathematical Modelling and Analysis 

Let us consider the peristaltic motion of blood on a porous channel. Blood is treated as a 
viscous couple stress fluid (non-Newtonian). We take (X,Y) as Cartesian coordinates of a point, 
X being measured in the direction of wave propagation and Y in the direction normal to the 
mean position of the corresponding small blood vessel. In the analysis that follows we consider 
channel flow. The constitutive equations and equations of motion for couple stress fluid flow in 
the absence of body moment and body couple can be put as [41, 36, 37, 38] 

dvi 

r S ij = P-^ (1) 

e ijk Tf k + Mjij = (2) 
Tij = -PSij + 2/idij (3) 
Hij = 4r)u jti + Ar/'uij (4) 
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in which p is the fluid density, and designate the symmetrical and antisymmetrical parts 
of the stress tensor Ty, Vi the velocity vector, My the couple stress tensor, p^ designate the 
deviatoric part of My, stands for the vorticity vector, dij is the symmetric part of the velocity 
gradient, rj and rj are constants associated with the couple stress effect, P the pressure. 

Let Y = H and Y = —H be respectively the upper and lower boundaries of the channel 
(cf. Fig. 1). The motion is considered to be induced by a progressive sinusoidal wave train 
propagating with a constant speed c along the channel wall, so that we have 

H = d + acos(^(X -ct)), 
where a is the wave amplitude, A the wave length and d the half width of the channel. 
We shall use Brinkman model for investigating the flow through the porous medium. Thus the 
equations that govern the flow of the couple stress fluid in the porous channel may be written 
in the form 



in which k is the permeability parameter, 77 the couple stress parameter, e stands for the porosity 
parameter and V 2 = + • 

Considering a wave frame (x, y) that moves with a velocity c away from the fixed frame (X, Y), 
let us use the transformations 

x = X-ct, y = Y, u = U - c, v = V, p(x, y) = P(X, Y, t) (8) 

in which (u, v) and (U, V) are the velocity components in the wave frame and the fixed frame 
respectively, p and P stand for pressure in wave and fixed frames of reference. 
Therefore the equation governing the flow of the fluid which is steady in the wave frame of 
reference can be written in the form 

4rt)-| + r v2 *-' v4 *-^-f (10) 
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Let us now introduce the following non-dimensional variables: 

x _ y _ u v d d 2 p - ct H a 

A a c co A [icX Ada 

2 _V/V _ I 2 R _pcd k 

Here I 2 is a material constant that has dimension of length square. It can be identified with a 

property that depends on the size of the fluid molecule (cf. Singh [59]). Dropping the bars over 

the symbols, the equations governing the steady flow of the fluid can be rewritten as 

j.du dv q 
dx dy 

RS ( * + = + I V \ - a^u - ±u - ^ (14) 

y ox ay ) ox e Da Da 

y ox Oy J oy e Da 

d 2 d 2 

in which V 2 = 5 2 — + — (16) 
ox z oy 1 

Boundary Conditions: 

Keeping in view the physical conditions of the problem, the boundary conditions in the fixed 
frame may be put mathematically as 

{%) U Y = 0, Uyyy = at Y = 

an 

(ii) U = , V = — , -(V xx - U YX )H X + (V XY -U YY ) = at Y = H 

ot 

(17) 

In the wave frame, these boundary conditions in terms of the non-dimensional variables will be 
transformed to 

(l) Uy = 0, Uyyy =0 Ot J/ = 

dh 

(ii) u = 0, v = — , - (S 4 v xx - 5 2 u yx )h x + (5 2 v xy - u yy ) = at y = h 

(18) 

Applying long wave length approximation (5 <C 1), for small Reynolds number (cf. Shapiro et 
al [52]), the equations (14) and (15) reduce to 



and = (20) 
dy 



d 2 

where 



dy 2 

Thus the pressure turns out to be a function of x only. The solution of equation (19) satisfying 
the boundary conditions(18) is found in the form 



Da m\ cosh(miy) dp Da m\ cosh(rri2y) dp dp 
(ml — m\)cosh{mih) dx (m 2 — m 2 )cosh(m 2 h) dx dx 



where 



, l-Jl- 4e 2 a 2 /Da „ 1 + Jl- 4e 2 a 2 /Da 

-h<y<h,m\ = V - 2 and m = V - 2 (22) 

2e 2e 



The rate of volume flow 'q' through each section is a constant (independent of both x and t). 
It is given by 

rh 

q= u dy. 
Jo 

Substituting (21) and performing the integration, we find 

Da m\ tanh[m\h) dp Da m\ tanh(m 2 h) dp ^ ^dp ^ 
mi(m 2 — m 2 ) dx m2{m^ — m 2 ) dx dx 

Hence the flux at any axial station in the fixed frame is found to be given by 

,-h 

Q= (u+ l)dy = q + h , 
Jo 

while the expression for the time-averaged volumetric flow rate over one period T(= of the 
peristaltic wave is obtained as 

Q = fJ*Q dt = f£( ( i + h ) dt = Q + 1 ■ 

The pressure gradient obtained from equation (23) can be expressed as 
dp (q + /j)m 1 m 2 (m2 — m\) 



dx Da m\ tanh[m\K) — Da m\ tanh(m 2 h) — Da hmim 2 (m2 — m 2 ) 
The pressure rise per wavelength and the frictional force can be calculated by using the relations 

Jo V dx ) 
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dx. 



Integrating over one wave length, we find 



mim2(m| — rrify ■= mim^m^ — mf) 

= f; Qh H f; ^2 

.Da Da 

m\m2(m\ — m\)- m\m2(m\ — mX) 
F = Qh H h 



Da 



Da 



where 



Jo 



da; 



/ 2 = 
/ 3 = - 

/o tanh(m\h) — m\ tanh(m 2 h) — /imimjfm^ — m\) 
In the fixed frame of reference, the expression for the non-dimensional axial velocity reads 



m\ tanh(mih) - 


m\ tanh(rri2h) - 


■ hmii7i2(m2 — mf) 


1 


(h - l)dx 




m\ tanh(m\h) 


— mf tanh(m,2h) 


- /im 1 m 2 (m2 — mf) 


1 


hdx 





m\ tanh(m\h) — m\ tanh(iri2h) — /jm 1 m 2 (m2 — m i) 
1 h(h - l)dx 



U(X,Y,t) = 



m 1 rri2(Q + h - 1) 



171% tanh(mih) — m\ tanh(m 2 h) — /im 1 m 2 (m2 — mf) 

( ml cosh{m\Y) m? cos/i(m 2 V) 9 9 \ 
' + m\ - m l 2 \ 



x 



cosh[m\K) cosh(m2h) 
with h = \ + <f) cos[2ti(X - t)] 




♦ 



— Study of Shapiro et. al. [25] 
• • • • Our study 

0.2 0.4 0.6 0.8 f 



Fig. 2.1: Variation of Q with for Ap=0, a = 0.001, Da=10000, e=0.9999. 
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3 Results and Discussion 



Analytical expressions of the pressure gradient, the volumetric flow rate, force of friction and 
axial velocity for the problem under consideration have been calculated and presented in the 
previous section. For the complexity of the problem, it was not possible to determine the 
analytical expression of Ap and F in terms of Q and other parameters. Therefore, we had 
to employ appropriate numerical methods along with the use of the software mathematica. 
Similarly, an analytical treatment was found inadequate for finding the velocity as a function of 
Ap, X, Y, t and other related parameters. This necessitated calculating first 

q DaAp — toiTO2(to| — m\)l2 
mim^ivn^ — mf)Ji 

For the computation of this also, we had to make use of numerical quadratures. In this 
section, the problem will be investigated further for a specific situation by employing appropriate 
numerical integration techniques with the help of mathematica software. For the quantitative 
study, the following ranges of values for different parameters involved in the model analysis have 
been considered : 

= 0.1 to 0.9; Ap = -40 to 40; a = 0.001 to 1.0; Da=0.05 to 10000, e=0.3 to 1.0, 
Q = -1 to 3 (Cf. Refs.[20, 56, 57, 58]). 

3.1 Pumping Characteristics 

Plots in Figs. 2 illustrate the variation of the volumetric flow rate of peristaltic waves with 
pressure gradient for different values of the amplitude ratio, couple stress parameter, perme- 
ability factor and porosity parameter. Fig. 2.1 shows that the results computed on the basis of 
our study in the absence of permeability, porosity and couple stress parameters perfectly match 
with the results reported by Shapiro et al [52] for Q max when Ap=0, Q = Q/ac. One may ob- 
serve from Fig. 2.2 that in the range of values of the pressure gradient examined in the present 
study, the volumetric flow rate increases with the increase in the amplitude ratio in the entire 
pumping region (Ap >0), in the free pumping zone (Ap =0) as well as in the co-pumping region 
(Ap <0) for Ap > —2.0. However, the trend reverses as soon as the pressure gradient drops 
below -2.0. These observations, of course, hold good for a particular set of values indicated in 
Fig. 2.2. Fig. 2.3 shows that in the entire pumping region the volumetric flow rate increases 
with the increase in couple stress parameter, whereas in the co-pumping region, a reverse trend 
is noticed. From the same figure, we further derive the information that for a Newtonian fluid 
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Fig. 2.4 Fig. 2.5 

Fig 2.2-2.5: Variation of Q with Ap (Fig 2.2) for different values of along with 
a=0.2, Da=10000, e=l; (Fig 2.3) for different values of a with 0=0.5, Da=10000, e=l; (Fig 
2.4) for different values of Da with 0=0.5, a=0.2, e=0.95; (Fig 2.5) for different values of e with 
0=0.5, a=0.2, Da=0.5 
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Fig. 3.1 (t 



= 0.0) 



Fig. 3.2 (t = 0.25) 




Fig. 3.3 (t = 0.5) Fig. 3.4 (t = 0.75) 

Figs. 3.1-3.4: Aerial view of the velocity distribution at different times 
(a = .003, Da = 10000, e = 1, Ap = 0, = .5) 
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(a ~ 0), the magnitude of the pressure rise is less than that for the couple stress fluid both in 
pumping and co-pumping regions. One can observe from Fig. 2.4 that the volumetric flow rate 
can be gradually reduced in the pumping region when Ap > 5 and enhanced in the pumping 
region when Ap < 5, free pumping region as well as in the co-pumping region by increasing 
the Darcy number through some mechanism. It is found from Fig. 2.5 that with a rise in the 
porosity parameter, volumetric flow rate increases in the pumping region, but reduces in the 
co-pumping region. 

3.2 Velocity Change 

An idea of the distribution of axial velocity can be from Figs 3 in the cases of free pumping (Figs. 
3.1-3.9), co-pumping (Fig. 3.10 and 3.12) and pumping (Fig. 3.11) zones for different values of 
the amplitude ratio, the couple stress parameter, the Darcy number and the porosity parameter. 
In the fixed frame of reference, since the velocity profiles and the height of the channel change 
with time, the axial velocity has been studied at the particular time instant t=T/4. Aerial view 
of the velocity distribution of a Newtonian fluid on the basis of our study has been presented in 
Figs 3.1-3.4, while Fig. 3.5 gives the velocity distribution of a couple stress fluid (a = 0.4) in 
the plane of the channel. We observe that at any instant of time there exists a retrograde flow 
region. But the forward flow region is predominant here as the time averaged flow rate Q is 
positive. Our study also reveals that for a couple stress fluid, there exist two stagnation points; 
for example, at t=0, one of them lies between X=0.25 and X=0.5, and the other between X=0.5 
and X=0.75. A similar observation was reported earlier for a similar study [53]. Figs 3.6-3.7 
reveal that for free pumping case (Ap = 0), the magnitude of the velocity for Newtonian as 
well as couple stress fluid in the forward and retrograde flow regions increases when the value of 
the amplitude ratio increases. Moreover, from Fig 3.7, one may observe that as the value of a 
increases, the velocity in the forward flow region increases. It may also be noted that the couple 
stress parameter a enhances the flow velocity in the retrograde flow region. The influence of 
the Darcy number Da on the velocity profile for a couple stress fluid (with a — 0.1) is depicted 
in Fig. 3.8. One can further note that as the Darcy number Da decreases, the axial velocity 
diminishes in forward region, whereas it increases in magnitude in retrograde flow region (i.e 
retrograde flow enhances). It is worthwhile to mention here that the parabolic nature of the 
velocity profile is disturbed in the central region for small values of Da. We further observe 
that for a couple stress fluid, the axial velocity decreases with an increase in the value of the 
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Fig. 3.7 
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Figs. 3.5-3.12: Longitudinal velocity profiles in different situations 
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porosity parameter e of the channel (cf. Fig3.9) in both forward and backward flow regions. 
Figs. 3.10-3.11 give the distribution of velocity of a couple stress fluid for co-pumping (Ap < 0) 
and pumping (Ap > 0) cases. From both the figures, it is revealed that as the value of Ap 
decreases, the velocity increases in both the flow regions. It is important to mention that unlike 
free pumping case, in the co-pumping case with Ap = — 1, the velocity reduces with the increase 
in the value of the couple stress parameter a (cf. Fig. 3.12). 

3.3 Frictional Force 

The average rise in frictional force in the case of a couple stress fluid are calculated over one 
wave period. Figs. 4. (a) depict the variation of frictional force with time for different values 
of Q when <fi = 0.5. It can be observed that the effect of increasing the flow rate Q is to 
enhance the frictional force F. As expected, from Figs. 4(b), we find that F decreases with rise 
in amplitude ratio. Figs. 4(c) reveal that the magnitude of the frictional force for a Newtonian 
fluid is less than that in the case of a couple stress fluid and that as a increases, the magnitude 
of F enhances. The quantum of influence of the Darcy number on frictional force is shown in 
Figs. 4(d). These figures indicate that the magnitude of F diminishes when the value of Da 
reduces, while Figs. 4(e) show that with an increase in the porosity parameter e, the magnitude 
of the frictional force decreases. 

3.4 Streamline Patterns and Trapping 

Under certain conditions the streamlines on the center line in the wave frame of reference are 
found to split in order to enclose a bolus of fluid particles circulating along closed streamlines. 
This phenomenon is referred to as trapping, which is a characteristic of peristaltic motion. Since 
this bolus appears to be trapped by the wave, the bolus moves with the same speed as that of the 
wave. Figs. 5.1 illustrate the streamline patterns and trapping for different values of the couple 
stress parameter a corresponding to Q = 1.3, = 0.3, Da = 10000, e — 1. It is important to 
observe that the size of the trapped bolus increases as the couple stress parameter a increases. 
Formation of trapping zone and streamline patterns is depicted in Figs. 5.2 for different values 
of Darcy number Da, when a = 0.003, e = 0.8, = 0.5, Q = 1.4. These figures reveal that with 
a reduction in Da, the trapping zone decreases and it disappears completely when Da attains 
the value 0.08. Plots showing the effect of the porosity parameter e on trapping are presented in 
Figs. 5.3. One can observe that size of the trapped bolus decreases with increase in e and vanish 
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Figs. 4: Frictional force versus time for different cases 

(a) a = 0.2, = 0.5, Da = 10000, e = 1 

(b) a = 0.2, Q = 0.1, Da = 10000, e = 1 

(c) = 0.5, Q = 0.1, Da = 10000, e = 1 

(d) a = 0.2, = 0.5, Q = 0.1, e = 0.95 

(e) a = 0.2, = 0.5, Q = 0.1, Da = 0.5 
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Figs. 5.2: Streamline patterns for different values of Da with a = 0.003, Q = 1.4, 
0.8 



0.5, e 



when e=l. Figs. 5.4 indicate that size of trapped bolus increases as the flow rate increases and 
that with increasing flow rate, the bolus gets shifted towards the boundary more and more. 



4 Summary and Conclusion 

Peristaltic transport of a couple stress fluid in a porous medium has been the concern of the 
present investigation. The study particularly pertains to a situation when the Reynolds number 
is low and curvature of the channel is quite small. The primary motivation behind the study 
has been to have an insight into the physiological problem of peristaltic transport of blood in a 
segment of the micro-circulatory system, when the segment is in a pathological state, whereby 
the lumen of the particular segment turns into a porous structure. Emphasis has been paid 
to investigate the peristaltic pumping characteristics, blood velocity, frictional force, streamline 
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Figs. 5.4: Streamline patterns for different values of Q with a = 0.003, <fi = 0.4, Da = 
10000, e = 1 
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patterns and trapping as a function of the peristaltic wave amplitude as well as Darcy per- 
meability and porosity parameters. Although the analysis starts with an analytical procedure, 
suitable numerical techniques have been employed whenever analytical methods have been found 
inadequate. The results and observations of this study (reduced to some particular cases) are in 
good agreement with those reported earlier by other investigators (cf. [52]). A very significant 
conclusion that can be made on the basis of this study is that it is possible to enhance the 
pressure rise Ap as well as the peristaltic pumping performance by increasing the amplitude 
ratio/couple stress parameter and/or by reducing the porosity. The study further reveals the 
important fact that change in velocity is strongly affected by the wave amplitude, the couple 
stress effect and the porosity factor. The other important conclusion that can be drawn out 
of this study is that while flow reversal is enhanced due to increase in pressure rise in the free 
pumping case, it is possible to bring about reduction in flow reversal by increasing the couple 
stress effect /permeability. Moreover, the occurrence of trapping can be reduced/eliminated by 
decreasing the permeability of the medium through the application of some appropriate mecha- 
nism. On the basis of this study we can make an important conjecture that when the pressure 
rise Ap is negative, the couple stress effect creates resistance to flow, whereby the velocity is 
reduced; however, the flow is enhanced due to the couple stress effect in the case of free pumping 
Ap = 

The present study reveals that it is possible to increase both pumping and pressure by in- 
creasing the amplitude ratio and couple stress parameter and also by reducing the permeability 
(Darcy number). The present study reveals further that the velocity change is strongly depen- 
dent on 0, a and porosity. Lastly, the magnitudes of Ap, and a affect significantly the nature 
of flow reversal. 
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